High Harmonic Generation in SF 6 : Raman-excited Vibrational Quantum Beats 
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In a recent experiment (N. Wagner et al. on SF6, a high- harmonic generating laser pulse 
is preceded by a pump pulse which stimulates Raman-active modes in the molecule. Varying the 
time delay between the two pulses modulates high harmonic intensity, with frequencies equal to the 
vibration frequencies of the Raman-active modes. We propose an explanation of this modulation as 
a quantum interference between competing pathways that occur via adjacent vibrational states of 
the molecule. The Raman and high harmonic processes act as beamsplitters, producing vibrational 
quantum beats among the Raman-active vibrational modes that are excited by the first pulse. We 
introduce a rigorous treatment of the electron-ion recombination process and the effect of the ionic 
Coulomb field in the electron propagation outside the molecule, improving over the widely-used 
three-step model. 

PACS numbers: 



High harmonic generation (HHG) is commonly under- 
stood as a 3 step process @] in which an electron ion- 
izes from a molecule, propagates in a strong laser field, 
and then recombines with the parent ion while emitting a 
photon. Acceleration by the laser field allows the electron 
to return with a large kinetic energy and emit photons 
with energy much higher than those of the driving laser. 

Although the HHG process is primarily electronic in 
character, recent experiments have shown that vibra- 
tional degrees of freedom can play a role. In the ex- 
periment by Wagner et al., [l( a 25 fs, 2.4 x 10 14 
W/cm 2 HHG laser pulse was preceded by a weaker 25 fs, 
5 x 10 13 W/cm 2 pulse which excited Raman-active vibra- 
tions in the molecule. The intensity of the high harmonic 
light was found to oscillate with the interpulse delay time 
at the excited molecular vibration frequencies. The sur- 
prising result was that the breathing mode, overwhelm- 
ingly dominant in conventional Raman experiments with 
this molecule, is no longer the strongest mode seen in the 
HHG experiment. 

In this Letter, we interpret these oscillations as an in- 
terference between indistinguishable quantum pathways 
associated with different intermediate vibrational states 
during the HHG process, illustrated in Fig. [T] We de- 
velop a quantum mechanical description of the recent 
pump-probe SFg experiments, using a framework that 
substantially improves on the three-step model. Q We 
have included the effects of the ion's Coulomb potential 
on the propagating electron, and utilized a nonpertur- 
bative electron-molecule scattering wavefunction [3j to 
calculate the recombination amplitude. The resulting 
calculation exhibits partial agreement with experimental 
observations. 

We adopt a level of approximation in which all op- 
erators depending on the nuclear coordinates are ex- 
panded to first order in the normal mode coordinates 
Qi, then the full vibrational state vector sepa- 



rates into a product of uncoupled normal mode vectors 
\yb{t)) [i) = af {t) |0) W + af (t) |1) W that can be treated 
individually. 

After a Raman pulse with the intensity and duration 
used in Ref. 1], only the |0) and |1) states of a given 
normal mode have significant amplitude. The following 
two-state, one-dimensional picture shows how individ- 
ual Raman-active modes affect the high harmonic signal. 
Atomic units are used throughout this work. 

During the Raman pulse, stimulated Raman scattering 
changes the vibrational state of the molecule from an 
initial |0) into a coherent superposition a |0) + ai |1) of 
the zeroth and first vibrational states. The vibrational 
coefficients follow equations of motion given by 

ia ni (t) =uji (m + -) a Hi (t) - ^^2E A (t)E B (t)x 

^ ' A,B (1) 

[aABd ni + diOtAB^Jrn + la„ i+ i + \JrHa ni -i)\ • 

Here u>j is the normal mode frequency, indices A and B 
run over {a:, y, z), EA(t) is the component of the electric 
field in the (body-frame) A direction at time t, Qi is the 
normalized displacement associated with normal mode 
i and cm_b(Qi, Q2, ■■■) is the polarizability tensor of the 
molecule. These equations of motion have off-diagonal 
elements only if diCtAB = (2mWj) dctAB /dQi\o 7^ 0, 
which is the condition for a mode to be Raman active. 
The polarizability tensor and its derivatives are found b 
performing an unrestricted Hartree-Fock calculation 
using the aug-cc-pVTZ basis set. [5[ 

Between laser pulses, |0) and |1) evolve as simple har- 
monic oscillator eigenstates, becoming i/vib = ao |0) + 
a\e~ iuJT |1) at the beginning of the high harmonic pulse 
for an interpulse delay of r, where the normal mode index 
(i) is omitted for brevity. Before the electron tunnels free 
of the ion, the high harmonic pulse stimulates the normal 
mode further according to Eq. |T]) . This is approximated 
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by a unitary 2x2 transfer matrix M_, where M nm is the 
amplitude to end in state \i) after starting in state \j) at 
the beginning of the pulse. 

Ionization and recombination, both electronic pro- 
cesses, are both modulated strongly by nuclear mo- 
tion. Taylor-expanding the tunncl-ionization operator / 
to first order in Qi, with the substitutions IqI = 7| eq 

and = (2muji) di/dQi\ cq and the identity Qi = 
(2muji) 1 ' 2 (Ai+A\), yields the first-order expansion into 
raising and lowering operators, / = lot + T,il[ l \Ai + 
A\). Here I is the identity operator, and ^4j,^4j are the 
raising and lowering operators for the i— th normal mode. 
For each mode considered, the reduced mass m is equal 
to the mass of a single fluorine atom. The recombination 

operator R can be derived to first order using identical 
logic. In both cases, dependence on nuclear positions 
means that the vibrational state changes along with the 
electronic state. 

The evolution of the ionic vibrational/electronic wave- 
function between ionization and recombination is in gen- 
eral quite complicated, since the three degenerate orbitals 
of SF^, which have Ti g symmetry, are coupled by vi- 
brational degrees of freedom [1, 0, S|- The linear and 
quadratic terms in the Jahn- Teller Hamiltonian @, 0, S] , 
which governs the coupled vibrational/electronic evolu- 
tion, are found for each Qi by fitting the eigenvalues of 
the coupling matrix to the lowest 3 adiabatic energies 
of SFji" for different displacements of the ion away from 
the maximum symmetry configuration. The energies are 
found using Gaussian's CASSCF method and a cc-PVDZ 
basis set [4]. In the notation of [6], Vr 2g =.001209 H/bohr, 
Ve 9 =.1420 H/bohr, iVi=-.0362 H/bohr 2 , P T2s = 7288 
H/bohr 2 , P Ss = 1.8486 H/bohr 2 . For the A lg mode, 
which does not enter into the vibronic Hamiltonian, an 
adiabatic potential E — VA lg QA lg + l/2^Ai g Q 2 i 1 > w ^h 
^=.0645 H/bohr, 70^=2.98 H/bohr 2 gives the po- 
tential energy surface for all three electronic states. An 
important simplification is that the off-diagonal coupling 
between different electronic states, proportional to Vt, is 
small and can be neglected for the short times between 
ionization and recombination. 

Neglecting off-diagonal coupling between electronic 
states, the adiabatic potential felt by the ion in a par- 
ticular electronic state is Hi = J^- + ViQi + ^KiQ 2 , 
which can be rewritten to first order in the basis of 
oscillator states of the neutral molecule. Evolution of 
the vibrational state is given by a transfer matrix TV = 
exp(-iH(t rct - t ion )). 

In a two-state treatment, the i— th vibrational wave- 
function of the neutral molecule after recombination has 
occurred is \tp V ib) — d |0) + d\ |1), where 

(do di) = (ao a ie ~ luJT ) M T / T A T if (2) 

and .A T is the transpose of matrix A. 



The number of photons emitted in a given harmonic 
is proportional to do ■ dg + d\ ■ d*. The high harmonic 
intensity is a sum over all Raman active modes i: 

P(t) = P Q + EjP^ cos (w iT + St) (3) 

. The static Pq primarily results from terms of the form 
a^ao, while Pi results from terms of the form aoa*e lWT 

and a ia ^ iuJT . Defining W = P if - RNIM_ , P = 
agWooao and Pi cos (wt + <5) = ±(ale luJT W 10 a + c.c). 
Since I\ and Pi are small relative to I and R , only 
their first-order terms are kept. 

At this level of approximation, calculating / and R as 
functions of the nuclear coordinates and substituting the 
expectation values for Qi at ionization and recombina- 
tion would give identical results. Nevertheless, tracking 
the quantum mechanical pathways in this manner is in- 
formative, because it allows the prediction of other ob- 
servables less amenable to a "classical nuclear motion" 
analysis, like the relative populations in |0) or |1) after 
recombination. 

This framework is applied to real molecules, using an 
improved version of the 3 step model. For the ionization 
step, a simple one-dimensional WKB tunneling picture 
describes an electron tunneling only in directions parallel 
to the laser electric field. This is motivated by the "ini- 
tial value representation" [^, [l(| ■ In the classically for- 
bidden region under the barrier formed by the molecular 
potential and the laser field, the tunneling wavefunction 
equals the value of the unperturbed molecular HOMO at 
the inner turning point times a declining WKB exponen- 
tial. With the direction of the electric field as —z, the 
wavefunction at the outer and inner turning points are 
related by 

ipt(x,y,zt P 2,t) = ihiouo(x,y, z tp i,t)x 
|Ci/C 2 | 1/6 Bi(0)/Ai(0)ex P [- / dzk(z)], {l 

J Ztpl 

where k(z) — y/2(V mo \ + Masor — E), Ai and Bi are Airy 
functions, ztpi and zt P 2 are the inner and outer turning 
points, C\ and C2 are the z components of the slopes of 
2V{r, t) at the two turning points, and the path integral 
is calculated along the z direction parallel to the applied 
electric field. 

After tunneling, the free electron wavefunction's evo- 
lution is relatively simple until it rescatters from the par- 
ent ion. ip c {r, O, t), the continuum wavefunction at the 
instant just prior to the electron rescattering from the 
molecular i on, is found using Gutzwiller's semiclassical 
propagator [111 ]: 

K(r,t;r ,t Q ) = (2iri)- 3 / 2 ^/C(r,t;r Q ,t Q )x 
e-xp[iS(f,t;r ,t ) - i4>] 

Here S(r, t; rb, to) is the action integral S = J L(q, q, t)dt 
calculated for the classical trajectory starting at (r, t) and 
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is the 



ending at (f , t Q ) and C(r, i; f ,t ) = \- dro AdrE 
density of trajectories for given initial and final points. <j) 
is a phase factor equal to § times the number of conju- 
gate points crossed by the trajectory. The semiclassical 
continuum wavefunction is 

4> c {r,t) = / d 3 f / dt K(r 7 t;r ,to)ipt{ro,to) (6) 



When the electron recollides with the parent ion, its 
wavefunction is distorted strongly by the molecular po- 
tential and by exchange effects with the other electrons, 
which can dramatically change amplitudes for recombi- 
nation with respect to the plane wave approximation. 
Techniques described in Refs [3, 12, l]| determine a com- 
plete set of stationary field- free electron- molecule scatter- 
ing states. Beyond the range of the molecular potential, 
the scattering states are given in terms of incoming and 
outgoing Coulomb radial functions f Ei (r) and the scat- 
tering S-matrix as 



i>E,i m (r) =r-j=f El {r)Y l , r , 



t4 E StAr)y V ,m^^)Sl',m';l, m {E). 



(7) 



During the electron-molecule scattering, when re- 
combination occurs, the electron wavepacket is ex- 
panded in terms of these scattering states as ips = 
S dE Y,i, m A i,m( E )^E,i m e- lEt , where 



Ai, m (E) = e 



iEi 



d 3 ri> E , m (r)Mr,t). 



(8) 



For a chosen time of projection onto the scatter- 
ing states, Eqs. ((6]) and ([8]) together define a seven- 
dimensional integral over initial and final positions and 
initial times. However, the integrand oscillates rapidly 
almost everywhere, causing cancellations. Stationary 
phase techniques identify the region where the integrand 
oscillates slowly, which permits evaluation. 

The semiclassical action S is expanded to second or- 
der around the starting and ending points of some clas- 
sical trajectory, and about the starting time. Near 
the starting position and time, tpti^o + 8ro,to + Sto) = 
^t{r ,to)exp[i(k 5r - iEHOMO$t )}, while f E> i(r + Sr) = 

fEA r ) ex P[Tk E ,i(r)Sr], where k E ,i{r) = V / 2(E - V t (r)), 
near the final position. Angular derivatives of the spher- 
ical harmonics are neglected. 

The integrand oscillates most slowly when its complex 
phase is nearly constant. This happens when the linear 
terms in the 5* expansion are canceled by the linear terms 
in the complex phases of the initial and final states. For 
such a trajectory, the contributions of nearby trajectories 
with nearly the same XQ,t^,x,t to the integral will tend 
to add constructively. These "stationary phase" trajec- 
tories are not the only trajectories of interest in the prob- 
lem, but expanding S about their beginning and ending 



points describes the phase in the slowly- varying region of 
interest. 

These stationary phase conditions are met by a tra- 
jectory that begins "downstream" of the molecule in the 
direction of the electric field with zero momentum, moves 
only radially and parallel to the electric field until it 
reencounters the parent molecule at the scattering state 
energy. Only the incoming-wave part of the scattering 
states gives a nonvanishing contribution to the expansion 
coefficients. Performing the resulting gaussian integrals 
about the initial and final points, and about the initial 
time, yields expansion coefficients 



A l>m (E) 




dr A 



dpo,B 

i[S(r,t,r ,to) 



/^W* (*) 



(9) 



] Mr )e l ~ lEHOMoto] 



for a stationary phase trajectory that starts at (zo,£o) 
and ends at (z,t), where z is the direction of the electric 
field at the time of ionization. 

The recombination amplitude is now D(E) = 

J2i,m A i,m( E ) d Lm(E), where d E ,i m = (4>g\x\ip E ,i m } is 
the recombination amplitude calculated for each individ- 
ual scattering state. This connects to the quantum paths 
framework when / = ipt(ro, to), and R=D(E) calculated 
for J = 1. Both quantities are calculated at the equilib- 
rium geometry, then for a geometry distorted by 0.1 bohr 

in the normal mode coordinate to find I,R and their de- 
rivitives. This involves recalculating the scattering states 
and recombination dipoles for the distorted geometry. 

The modulation of the 39th harmonic, which the JILA 
experiment considered in detail, is calculated for com- 
parison with experiment. This harmonic is close to the 
measured cutoff, so it can only be produced by a laser 
half-cycle where the maximum of the electric field falls 
close to the maximum of the Gaussian envelope. Accord- 
ingly, only the modulation for a single half-cycle where 
these two maxima coincide was used. In the expression 
for the harmonic intensity at some harmonic order, Eq. 
([3]), all the quantities depend on molecular orientation. A 
rotational average of Pq and Pi were calculated, since the 
JILA experiment was performed on a gas jet of molecules 
which had no preferred orientation. Only polarizations 
perpendicular to the propagating beam were included in 
these averages. Although each of the T2 g and E g modes 
modulates the harmonic intensity strongly at particular 
orientations, the phase offset S in Picos(uit + 5) changes 
with orientation, cancelling some of the observed oscil- 
lation. The more symmetric A g mode experiences less 
cancellation because the initial Raman pulse stimulates 
it equally for all molecular orientations. The trends in 
the peak intensities of different normal modes are, how- 
ever, similar with and without the phase. Fig. [2] com- 
pares the spherically averaged peak-to-peak modulation 
with phase information included, with phase information 
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excluded by setting 5 — for all orientations, and the 
experimentally measured modulations for each mode for 
the two runs for which all three mode modulations could 



be distinguished from the background [14|. It should be 
noted that although the relative peak intensities are not 
in agreement with experiment, the intensities for the dif- 
ferent peaks are all the same order of magnitude, whereas 
in the Raman spectrum the A g mode is 20 times more 
intense than the others. [H 

Since the scattering wavefunction is expanded in field- 
free scattering states, the calculated modulation varies 
slightly depending on the time at which the semiclassi- 
cal wavefunction is projected onto the scattering states. 
Additional uncertainty may arise because of the various 
approximations made in the treatment of molecular scat- 
tering, detailed in 

B E3, El ■ Fig. [U shows modulations 
calculated when this projection is made at tot = 3.9 ra- 
dians, when the short electron trajectories return to the 
vicinity of the parent molecule with the correct energy to 
yield 39th harmonic photons upon recombination. 

In the present work, we have combined for the first 
time a rigorous treatment of electron-molecule scattering 
with a semiclassical description of electronic propagation, 
resulting in a flexible and robust implementation that 
has been used to treat a complex molecule with many 
internal degrees of freedom to a level of sophistication 
that is unprecedented in this area of research for such a 
large system. We have shown how the internal degrees of 
freedom allow for phenomena which have no analogue in 
atomic systems, demonstrating that the HHG signal can 
be modulated interferometrically by a molecule's vibra- 
tional state. Thus HHG could serve both as an interfero- 
metric probe of chemical dynamics (similar to the way in 
which rearrangement dynamics are investigated in Ref. 
[IH) and as a complementary tool to traditional spec- 
troscopic techniques, offering great potential for future 
investigations. 
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FIG. 1: (Color online) The vibrational interference model in 
one dimension. The molecule ends the Raman pulse in a su- 
perposition of the v=0 and v=l vibrational states. After a 
time delay, the two vibrational states are mixed by stimulated 
Raman scattering, "hopping" during ionization and recombi- 
nation, and evolution of the ionic wavefunction while the elec- 
tron is away. Interference between adjacent vibrational states 
modulates the hieh harmonic sienal. 
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FIG. 2: (Color online) Peak-to-peak modulation of the high 
harmonic signal vs. wavenumber, comparing theory to the 
two experimental runs for which data is available. Modula- 
tions corresponding to the same frequency have been placed 
side-to-side for purpose of comparison. 



(1985). 

[7] W. Moffitt and W. Thorson, Phys. Rev. 108, 1251 
(1957). 

[8] I. Bersuker, The Jahn-Teller effect and vibronic interac- 
tions in modern chemistry (Plenum Press, 1984). 
[9] W. H. Miller, J. Phys. Chem. A 105, 2942 (2001). 
[10] H. Nakamura, J. Theor. Comp. Chem. 4, 127 (2005). 
[11] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer, New York, 1990). 
[12] S. Tonzani and C. H. Greene, J. Chem. Phys. 124, 054312 
(2006). 

[13] S. Tonzani, Comp. Phys. Comm. 176, 146 (2007). 
[14] N. Wagner, private communication. 

[15] S. Baker, J. S. Robinson, C. A. Haworth, H. Teng, R. A. 
Smith, C. C. Chirila, M. Lein, J. W. G. Tisch, and J. P. 
Marangos, Science 312, 424 (2006). 



